Monodromies and functional determinants in the CFT driven quantum cosmology 
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We apply the monodromy method for the calculation of the functional determinant of a special 
second order differential operator F — —d 2 /dr 2 + g/g, g = d 2 g/dr 2 , subject to periodic boundary 
conditions with a periodic zero mode g — g(r). This operator arises in applications of the early 
Universe theory and, in particular, determines the one-loop statistical sum for the microcanonical 
ensemble in cosmology generated by a conformal field theory (CFT). This ensemble realizes the 
concept of cosmological initial conditions by generalizing the notion of the no-boundary wavefunction 
of the Universe to the level of a special quasi-thermal state which is dominated by instantons with 
an oscillating scale factor of their Euclidean Friedmann-Robertson- Walker metric. These oscillations 
result in the multi-node nature of the zero mode g(r) of F, which is gauged out from its reduced 
functional determinant by the method of the Faddeev- Popov gauge fixing procedure. The calculation 
is done for a general case of multiple nodes (roots) within the period range of the Euclidean time 
r, thus generalizing the previously known result for the single-node case of one oscillation of the 
cosmological scale factor. The functional determinant of F expresses in terms of the monodromy 
of its basis function, which is obtained in quadratures as a sum of contributions of time segments 
connecting neighboring pairs of the zero mode roots within the period range. 

PACS numbers: 04.60. Gw, 04.62,+v, 98. 80. Bp, 98.80.Qc 
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I. INTRODUCTION 

Successful calculation of quantum corrections on non- 
trivial backgrounds always presents a challenge and can 
be accomplished in a closed analytic form only in a lim- 
ited number of cases. This applies also to calculations 
in the quantum mechanical sector of field models even 
despite simplifications occurring in this sector due to its 
spatial homogeneity or other symmetries. A particular 
case of such calculations is the class of problems in cos- 
mology involving the differential operator of the form 



Id 2 d 1 
g dr ® dr g 



dr 2 
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(1) 



where g — g(r) is a rather generic function of the cosmic 
time depending on the behavior of the cosmological scale 
factor a and its time derivative, a = da/dr. From cal- 
culational viewpoint, the virtue of this operator is that 
<?(t) represents its explicit basis function - the solution 
of the homogeneous equation, 



Fg(r) = 0, 



(2) 



which immediately allows one to construct its second 
linearly independent solution and explicitly build the 
Green's function of F with appropriate boundary con- 
ditions. On the other hand, from physical viewpoint 
this operator is interesting because it describes long- 
wavelength perturbations in early Universe cosmology, 
including the formation of observable CMB spectra [Jill, 
particle creation in cosmology etc. In particular, 
for superhorizon cosmological perturbations of small mo- 
menta k 2 <§; g/g their evolution operator only slightly dif- 
fers from ((T|) by adding k 2 to its potential term, whereas 
in the minisuperspace sector of cosmology, corresponding 
to spatially constant variables, the operator has exactly 
the above form. 



Up to an overall sign, this operator is the same both in 
the Lorentzian and Euclidean signature spacetimes with 
the time variables related by the Wick rotation r = it. 
In the Euclidean case it plays a very important role in 
the calculation of the statistical sum for the microcanon- 
ical ensemble in cosmology. This ensemble realizes the 
concept of cosmological initial conditions by generalizing 
the notion of the no-boundary wavefunction of the Uni- 
verse [H to the level of a special quasi-thermal state and 
puts it on the basis of a consistent canonical quantization 
p, @ . This concept is very promising from the viewpoint 
of the cosmological constant, inflation and dark energy 
problems In particular, it provides a thermal input 

in the red tilt of the COBE part of the CMB spectrum 
[Io| additional or, perhaps, even alternative to the con- 
ventional one generated from primordial vacuum fluctu- 
ations in early inflationary Universe . 

In this statistical theory context, when r plays the role 
of the Euclidean time, the properties of this operator es- 
sentially differ from the Lorentzian dynamics. In the lat- 
ter case the function g is a monotonic function of time 
because of the monotonically growing cosmological scale 
factor, whereas in the Euclidean case g(r) is periodic just 
as the scale factor a(r) itself and, moreover, has zeroes 
at turning points of the Euclidean evolution with a = 0, 
because g(r) oc d(r). This does not lead to a singular 
behavior of F because g also vanishes at the zeroes of 
9 [H, Q ! and the potential term of (JXJ) remains analytic 
(both g and g simultaneously have a first-order zero as it 
follows from the dynamical equation for a(r) in Q). Nev- 
ertheless, the calculation of various quantities associated 
with this operator becomes cumbersome due to the roots 
of g(r). Among such quantities is the functional deter- 
minant of F which determines the one-loop contribution 
to the statistical sum of the CFT driven cosmology of 
Since this operator has an obvious zero mode which 
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is the function g(r) itself, the functional determinant of 
F should, of course, be understood as calculated on the 
subspace of its nonzero modes. Thus, the focus of this 
paper is the calculation of such a restricted functional 
determinant of ([T]), denoted below by Det*.F. This de- 
terminant is calculated on the space of functions periodic 
on a compactified range of the Euclidean time r (forming 
a circle) with its zero mode g removed or gauged out. 

There exist several different methods for restricted 
functional determinants. When the whole spectrum of 
the operator is known this is just the product of all non- 
zero eigenvalues. With the knowledge of only the zero 
mode, one can use the regularization technique [l2[ or 
contour integration method [HI, H3] to extract the reg- 
ulated zero-mode eigenvalue from the determinant and 
subsequently take the regularization off. Here we use an- 
other approach to the definition of Det*.F dictated by the 
auge-fixing procedure for the path integral in cosmology 
. This procedure looks as follows. 
As shown in [6j, for a spatially closed cosmology with 
5 3 -topology the microcanonical statistical sum is given 
by the path integral over the metric g^ v and matter 
fields which are periodic on the Euclidean signature 
spacetime with a compactified time r (S 1 x S 3 topol- 
ogy). This integral can be transformed by decompos- 
ing (g fll ,(x),<j)(x)) into the minisuperspace sector of the 
lapse function N(t) and the scale factor a(r) of the 
FRW metric and inhomogeneous "matter" fields <P(x) = 
(4>(x), ip(x), A^x), h^(x), ...), 

Z=J D[g^4>] e ~^,H = J D[ajN] e -r[«,iV] )( 3 ) 

-r[a,N] = f m ^ e -S[a,N;0{x)] _ ^ 



Here, r[a, N] is the Euclidean effective action of the 
fields <P (which include also the metric perturbations h^ v ) 
on the FRW background, ds 2 = N 2 (t) dr 2 +a 2 (r) d 2 n^ 3 \ 
S[a, N; <P(x) ] = S[g^,(f>] is the original action rewrit- 
ten in terms of this minisuperspace decomposition and 
D[a,N] is the canonical measure of integration over 
periodic histories (N(t) 7 o,(t)^, which incorporates the 
Faddeev-Popov gauge fixing for the one-dimensional dif- 
feomorphism invariance of r[a, N]. 

In the theory with a primordial cosmological constant 
and a large number of free (linear) fields conformally cou- 
pled to gravity - conformal field theory (CFT) - such 
effective action is dominated by the contribution of these 
fields and is exactly calculable with the aid of their local 
conformal anomaly jT^-[l3 • The structure of the result- 
ing action [f| is 

r[a,N] = I dTN£(a,a')+F(r}), r, = <j> dr—, (5) 



where 



da/Ndr and the integration runs over the 



period of r on the circle S 1 of S 1 x S 3 . Here the effective 
Lagrangian of its local part C(a, a') includes the clas- 
sical Einstein-Hilbert term and the polarization effects 



of quantum fields and their vacuum (Casimir) energy. A 
nonlocal part of the action F(rf) is the free energy of their 
quasi-equilibrium excitations with the effective tempera- 
ture given by the inverse of the conformal time rj. 

The operator |T]) arises in the saddle point approxi- 
mation for the path integral ([3]). It enters the quadratic 
part of the action §5§ in perturbations Sa(r) and SN(t) 
on the background of the classical solution (a(r), AT(r)) 
for this action. In terms of the new set of perturbation 
variables 




SN 
a 



6 a 
^2 ' 



da'da' 



(6) 
(7) 



the part of the action quadratic in ip{r), which is the 
canonically normalized perturbation of the cosmological 
scale factor [9(, reads as 




dT(p(T)Fip(T), 



(8) 
(9) 



along with the zero mode g(r), which is given as the 
explicit function of the background history a(r) and its 
derivative ct(r). 1 

This mode serves as a generator of the global gauge 
transformation of the gravitational variable ip, S e ip(r) oc 
g(r) e. This transformation is the residual symmetry of 
the action ([5]) which remains after gauge-fixing the local 
diffcomorphism invariance of the original action ([5]) in 
the gauge h(r) = imposed on the second of the vari- 
ables ((6]) . Therefore, this symmetry is also subject to the 
Faddeev-Popov gauge fixing procedure which consists of 
imposing another gauge xi'P) — and inserting in the 
path integral the relevant Faddeev-Popov factor. This 
gauge condition x(v) an( i * ne Faddeev-Popov ghost fac- 
tor Q/H5II can be chosen in the form 

x = jdTk( T ) v {T), (10) 



dTk(T)g(r), 



g\Y= <bd T g\T), 



(11) 
(12) 



where fc(r) is a gauge fixing function and the role of the 
generator in the gauge transformation of v ? ( T ), $ El p{ T ) — 
R(r) e, is played by the the zero mode R(r) — <?(t)/||<?|| 
which is normalized to unity with respect to L 2 inner 
product on a circle. 2 



1 Here and in what follows we we choose the gauge N = 1 for the 
background solution (a(r), iV(r)) , so that time-parametrization 
invariant derivative a' = a/N = a. 

2 In local gauge theories the Faddeev-Popov path integral is not 
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Thus, integration over if takes the form of the Gaussian 
functional integral with the delta-function type gauge, 

(Det.F)" 1 / 2 

= const x J Dtp5{x) exp ^— i <j> drip F(p^j , (13) 

and serves as the definition of the restricted functional 
determinant of F. This definition is in fact independent 
of the choice of gauge by the usual gauge independence 
mechanism for the Faddeev-Popov integral. In particu- 
lar, enforcing the gauge x = means that the field (p is 
functionally orthogonal to the gauge fixing function k(r) 
in the L 2 metric on a circle, and the above definition is 
independent of the choice of this gauge fixing function. 

The calculation of the restricted functional determi- 
nant of ([T]) as an explicit functional of g(r) has already 
been undertaken in |18f for a particular case of this func- 
tion g(r) having only two zeroes in the period T of the 
time variable. Here we generalize this single-node result 
to the multiple nodes case with an arbitrarily high even 
number 2k of roots of g{r) within its period range. 3 This 
is the case of the CFT driven cosmology whose statisti- 
cal sum at different values of the primordial cosmological 
constant is dominated by the countable set of instantons 
having k oscillations, k = 1, 2, of the scale factor dur- 
ing the Euclidean time period [a @| _ the so-called gar- 
lands which carry the multi-node zero mode <?(t) oc a. 
First we express this determinant in terms of the mon- 
odromy of the second basis function ip(r) of F, which is 
linearly independent of g(r). In this way we actually re- 
produce the known monodromy formula of McKane and 
Tarlie for the restricted functional determinant of the op- 
erator subject to periodic boundary conditions [13 - [l4j . 
Then we derive the answer for this monodromy as an ad- 
ditive sum of contributions of segments of the time vari- 
able, Tj_i < t < Ti, i = 1,2, ...,2k, separating various 
pairs of neighboring roots of g on. S 1 , gin) = 0. Each 
contribution is given by a closed integral expression in 
terms of g(r) on an underlying segment. Finally, we re- 
produce the previously obtained result for a single-node 
case of [l8| and accomplish the paper with conclusions. 



invariant under arbitrary rescalings of gauge generators R. Their 
normalization is always implicitly fixed by the requirement of 
locality and the unit coefficient of the time-derivative term in the 
gauge transformation of Lagrange multipliers, Re ~ 1 X e + 
(which follows from the canonical quantization underlying the 
Hamiltonian version of the Faddeev-Popov path integral). For 
the residual global symmetry we do not have a counterpart in 
canonical formalism, and such a founding principle as canonical 
quantization does not seem to be available. Therefore, we choose 
this normalization with respect to L 2 unit norm corresponding 
to the canonical normalization of the variable tp in ((SJl . From the 
viewpoint of the definition of Det„.F as the product of operator 
eigenvalues, this corresponds to the omission of a zero eigenvalue 
of F. 

3 Since a periodic function can have only an even number of roots 
within its period, the case of the lowest nonvanishing number of 
roots we will call a single-node one. 



II. SETTING THE PROBLEM AND SUMMARY 
OF RESULTS 

Let us start with the setting of the problem and a 
brief formulation of the final result. The operator ([I]) is 
defined on a circle range of the time variable S 1 having 
the circumference T, which is parameterized by r 



r < r < r + T, 



(14) 



with the points tq and tq + T being identified. This range 
can be infinitely extended to the whole axis — oo < r < 
oo, multiple covering of S , on which the function g(r) 
and, consequently, the operator F are periodic with the 
period T, 



1 (r + T) = g(T). 



(15) 



In the fundamental domain this periodic function has 2k 
simple roots 



r < n < r 2 < ... < r 2fc = t 
9 {n) = 0, gin) + 0, 
g(n) = o, 



T, 



(16) 
(17) 
(18) 



one of them coinciding with the final (or starting) point 
of this domain. Another important assumption is that 
the second order derivative of this function at its roots is 
vanishing, which will be important for analyticity prop- 
erties of our formalism. 

For any two functions ipi and if2 the operator F deter- 
mines their Wronskian W[(p\, ^2] and Wronskian relation 
defined by the equations 



J dnp{F*(p 2 = J dT(p\F~y>2 — W[ipi,<p2[ 



T+ 



(19) 
(20) 



Arrows here denote the direction of action of the opera- 
tor F, i. e. ipjF — (F(pi), and the Wronskians appear as 
total derivative terms generated by integration by parts 
of the derivatives in F. When both ipi and ip2 satisfy a 
homogeneous equation with the operator F , their Wron- 
skian turns out to be constant. Also the vanishing Wron- 
skian implies linear dependence of these solutions. 

Let us now consider the solution of the homogeneous 
equation ip(r) normalized by a unit Wronskian with g 



F^{t) = 0, 
W[g,^] = l 



(21) 
(22) 



Together with <?(t) this solution forms a set of linearly in- 
dependent basis functions of F . However, in contrast to 
g{r) the basis function is not periodic, because we assume 
that the operator ([T]) has only one periodic zero mode 
smoothly defined on a circle (|T4")) . On the other hand, 
when considered on the full axis of r, due to periodicity 



4 



of g(r) this operator is also periodic F(t + T) = F{t). 
Therefore ip(r + T) is also a solution of the equation 
F(t)iP(t + T) = 0, and consequently it can be decom- 
posed into a linear combination of the original two basis 
functions with constant coefficients 



^(t + T) = V(t)+A 5 (t). 



(23) 



The unit coefficient in the first term follows from the con- 
servation in time of the Wronskian of any two solutions 
of the equation (|2ip . periodicity of g(r) and an obvious 
fact that W[g,g] = 0, 1 = W[g(r + T),^(t + T)\ = 
W{g(r), -0(r) + Ag(r)]. The coefficient A in the second 
term is nontrivial - it parameterizes the monodromy of 
-0(r) and will play a central role in the construction of 
the determinant. 

Below we show that in terms of this monodromy pa- 
rameter and the square of the norm of the zero mode 
(fT2j) the restricted functional determinant of the opera- 
tor F subject to periodic boundary conditions has a very 
simple form 



Det, F 



A <j> dr g 2 (r). 



(24) 



In fact, this is the McKane-Tarlie formula (Eq.(5.2) of 
fl2l |) obtained by the regularization and contour integra- 
tion methods [12Hl4l | based on the earlier work of For- 
man j!9| for a generic second order differential operator. 
We reproduce this formula by the variational method for 
functional determinants. Beyond this, the structure of 
the operator (fT]) makes the problem exactly solvable and 
allows one to find the monodromy (|23j) in quadratures as 
an explicit functional of g(r). The answer reads in terms 
of the set of solutions of (f2"Tj) explicitly defined by the 
following integrals for i = 1, 2k, 



&i(T)=g(T) 



dy 
5 2 (y)' 



r,-_i < r, t,- < 7V. 



(25) 



r where the integral (1251) is convergent because the roots 
of g{r) do not occur in the integration range. Its limits 
are well defined also at the boundaries of this segment, 



I 



g{n) 



(28) 



because the factor g{r) tending to zero compensates 
for the divergence of the integral at t — > Tj — and 
t — > Ti-x + 0. Moreover, because of (jT8j) the functions 
^(t) are differentiable in these limits, and all the quan- 
tities which enter the algorithm ([2"T]) are well defined. In 
particular, for any such time segment [Tj_i, t^] = [r_, t + ] 
the derivatives of ^(t) at its boundaries are given by the 
convergent integral (we omit the label i marking any of 
the functions 3^(t) and their auxiliary points t* in (f2"T|) ) 



*{t±) = / dy 



g{ T±)-g{y) 
g 2 (y) 



7(r*) 



(29) 



Note that the integrand here is finite at y —> t± because 
of g{j±) = 0. These properties of •Pi(r) guarantee that 
the obtained result is independent of the choice of the 
auxiliary point r* for each A^, and the monodromy pa- 
rameter l|26p is uniquely defined. 

Final comment concerns the overall normalization in 
flU]) . The formula of McKane-Tarlie [H-lli] or the vari- 
ational method which we use below, in principle, give 
only the ratio of determinants for two different operators 
(say with different functions g and periods T), whereas 
each determinant contains an infinite numerical factor 
generated by UV divergent product of eigenvalues of F. 
One can only say that this factor is independent of g(r) 
and T and depends on the UV regularization. Here we 
fix it by the zeta-function regularization in which it turns 
out to be equal to one. 



on the segments of r-range connecting the pairs of neigh- 
boring roots of g(r). Here r* are the auxiliary points 
arbitrarily chosen in the same segments, and all these so- 
lutions are normalized by the unit Wronskian with g(r), 
Wlg^i] = I. Then the monodromy parameter in ques- 
tion reads as the additive sum of contributions of all these 
segments 



2 k 



g{n) g{n-{) 



(26) 



(27) 



III. VARIATIONAL EXPRESSION FOR THE 
DETERMINANT 



Representing the delta function of the gauge condition 
in (|I3[) via the integral over the Lagrangian multiplier 
7r we get the Gaussian path integral over the periodic 
function <p(r) and the numerical variable tt, 



(Det*F) 



-1/2 



const x — - f D(pdn exp ( — S e s[<fi(T),w] 

\\g\\ J v 

const x Q DetFJ . (30) 



The main property of these functions ^(t) is that each 
of them is defined in the i-th segment of the full period of 



Here S c g [ <p(r), tt ] is the effective action of these variables 
and F is the matrix valued Hessian of this action with 
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respect to <P — (</?(r),7r), 

S cS [ip{t),tt} = <j> dr ^ i i^F(/3 - /,v/. 



IV. PERIODIC GREEN'S FUNCTION 



FS(t,t') -ifc(r) 
-ife(r') 



(31) 
(32) 



(note the position of time entries associated with the vari- 
ables (<^(t),7t) and = (<^(t'),7t)). 

The dependence of this determinant on g(r) and k{r) 
can be found from its functional variation with respect 
to these functions. From (l30l) we have 



8 In Det*F 



-28 In Q + 28 In | \g \ \ + Tr ( SF G ) , (33) 



where G is the Green's function of F, F G = I and the 
functional trace of any matrix with the block-structure 
of ([52")) is defined as 



Tr 



A(t,t>) B(t) 
B(t') a 



dr A(t, t) + a . 



The block structure of the matrix Green's function G 
has the form 



G(t,t') 

Q 







(34) 



where the Green's function G(r, r') in the diagonal block 
satisfies the system of equations 



FG(t,t') = S(t,t')- 
dT k(r ) G(t, t') =0 



Q ' 



(35) 
(36) 



which uniquely fix it. The second equation imposes the 
needed gauge, whereas the right hand side of the first 
equation implies that G(r, r') is the inverse of the oper- 
ator F on the subspace orthogonal to its zero mode. 

The trace of the functional block-structure matrix in 
(|33[) corresponding to the variation of g(r) reads 



Tr 



(<S fl FG) = Tr (SF G(t, r')j = I dr SF G(t,t') 



A similar variation of the gauge-fixing function gives a 
vanishing answer 

8k In ^Det» F^j = — ^ <j> dr 8k g 



-Tr 



lg{r)5k{T') 
W 1 
i§ dy8k{y)G{y 1 T l ) —§ dy 8k(y) g(y) 



0, (38) 



as, of course, it should be in view of the gauge indepen- 
dent nature of the Faddeev-Popov path integral. This 
guarantees the uniqueness of the definition of the reduced 
determinant Det* F. 



For the calculation of the variation (1371) above we 
need the Green's function of the problem (f3"5| - ([3"6"| which 
should be periodic on the circle (|14l) . This means that 
the value of G(r, r') and its derivative with respect to r 
should match at r = tq + and t = tq + T — 0. To achieve 
this property we will slightly extend the circle domain to 
the left of the point tq 



to 



<T <T +T, £ > 0, 



(39) 



with an arbitrarily small positive e and demand that the 
monodromy of G(r, r') should be vanishing for this small 
e-range of r near tq 

A t G(t,t') = G{t+T,t')-G(t,t') = 0, T -e < r < t . 

(40) 

The ansatz for G(t, t') can be as usual built with the 
aid of two linearly independent basis functions of the op- 
erator. One basis function coincides with the periodic 
zero mode g(r) and another one ip(r) can be composed 
of the set of functions ^(r) defined by (|25|) on various 
segments Tj_i < r < Ti. 4 For an arbitrary choice of 
auxiliary points r* in (|25[) the composite function 



tI>(t) = !^(t), 



r,:_i < r < t, 



(41) 



will be continuous in view of (|2"5|) , but the continuity of 
its derivative will generally be broken, because generally 
the equality <fj(r,— 0) = !?j+i(r,-+0) is not satisfied. How- 
ever, this equality for i = 1, 2, 2fe — 1 can be enforced 
by a special choice of these auxiliary points t* , becoming 
the equation for their determination. The solution for t* 
is unique, always exists and belongs to the corresponding 
segment Tj_i < t* < r^. 5 On the other hand, the conti- 
nuity of the derivative of iP(t) cannot be attained at all 
roots of <?(t), i = 1, 2, ...2k, because it would correspond 
to the existence of the second zero mode periodic on the 
circle, which is ruled out by construction. Therefore, the 
second basis function of F is not periodic on the circle, 
but in view of periodicity of the operator it satisfies the 
fundamental monodromy property (|23| . Below we con- 
struct the periodic Green's function of F in terms of the 
monodromy parameter A. 

The Green's function of the problem (|35 |) -([36 |) can be 
represented as a sum of the particular solution of the 
inhomogeneous equation (|35p and the bilinear combina- 
tion of the basis functions g and ip with the coefficients 



For the extended range {39) the missing ^o(t) can be denned by 
identifying r_i with T2k-i — T and choosing some Tg in T—\ < 
^ t * < r . 

5 Indeed, the quantity d&ifa) / dr* = —g(Ti)/g 2 (r*) is a sign defi- 
nite function of r* nowhere vanishing on this segment, its abso- 
lute value quadratically divergent to oo at its boundaries. This 
in its turn means that <Pi (r; ) is a monotonic function of r* which 
also ranges between — oo and +oo and therefore guarantees the 
unique solution for r* on this segment. 
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providing the periodicity property (1401) . In view of the 
Hermiticity of the operator F the kernel of its Green's is 
symmetric, G(r, r') = G(t',t), so that we will look for 
it in the form 



G(r,r') = G f (t,t') + - n(T,T')+aH^{Ty) 



(42) 



where 



G F (r, t') ee -{g{r)i,{r') ~ ^r)g(r')) 9{r r') 

+ k^(r)g(r') 9 (t)^{t')) 6(t' - r) , (43) 



SI(t,t') ee uj(T)g(T')+g(r)u(r') , 

H^(T,r') EE V(r)V(r') , 

H^{t,t') ee V(T)ff(r')+5(T)V(r') , 
# m (t,t') ee . 9 (r) 5 (r'), 



(44) 
(45) 
(46) 
(47) 



and the function w(r) will be defined below. The 
first term of (|4"2")l generates the delta-function in the 
right hand side of the equation (I55|) . the second 
term Vl{t,t')/Q gives -k(r) g(r')/Q, while H^(t,t'), 
Hipg ( r ; T ') an d Hgg( T > T ') represent symmetric solutions 
of the homogeneous equation with coefficients a, (3 and 
7 to be fixed by the periodicity condition and the gauge 
condition (l36|) . The concrete form of these functions can 
be explained as follows. 

The simplest form of a particular solution of 
FGf(t,t') = 5(t,t') can be written down with the 
aid of the Heaviside step function as the combination 
3 , (t)'0(t') 9(t — t') + (t <-> t') bilinear in two basis func- 
tions g and ip (normalized by the condition W[g, ip] = 1). 
But we will be interested in the kernel vanishing at the 
coincident arguments, which is easily attained by adding 
a special solution of the homogeneous equation, 

G f (t, t') = 9 {t)^{t') 6(t - t') + ^(r)g(r') 6{t' - r) 

-i(5(r)V(r')+^(T)ff(r')). 

This leads to (g3J). 

A particular solution of FSI(t,t') = — fc(r) g(r') has 
the form (|44[) where the function w(r) satisfies the equa- 
tion Flu(t) — —k(r). It has as the first integral 



W[gM{T) = Q{r) , 
Q(r) ee / dyg(y)k(y), 



(48) 



with an arbitrary r* . This can be further integrated with 
a particular boundary condition w(t„) = to give the 
final expression for w(r), 



u(t) = i]}(t) J dy g(y)k(y) - g(r) / dyip(y)k(y) 

= V>(r) Q(t) - 9 {t) I dy V>(y) k(y) . (49) 



Similarly to - 0( r ) this function has a nontrivial mon- 
odromy which can be derived from the monodromy of 

V>(r), 



t+T 



A T w(r) = ip(r+T) / dy g(y) k(y) - ip(r) dyg{y)k{y) 



T+T T \ 

,ir) | J - J) dy^(y)k(y) 

T 

^(r+T)Q + A. 9 (r) J dygk 

r* 

~5( T ) J dy^pk ~ g{r) j dy (i/j + Ag^jk 

T, T, 

r 

+g(r) J dyi\)k 

T» 

,P(T + T)Q-g(T) fdy^(y)k(y). (50) 



Thus finally 



A t ui{t) =V(T)g+g(r) -/ dy^(y)k(y) + AQ 



r,+T 



Using this monodromy and (|23[) we have 



(51) 



A r G F (r, 7-0 = i( ,g(r + I>(rO - V(t + T)ff(rO ) 
~(i>(r)g(r') - g(r)i;(r')) 
= g{r) i,{r') - </>(r) <?(t0 - ± A g{r) g(r') , (52) 

where we took into account that for the range of r in the 
small e-domain of the enlarged range (|39[) . tq— e < t < To 
and ro < t' < r + T the arguments of G f (t 7 t') satisfy 
the relations r < r' and r + T > r' (the latter inequality 
can always be enforced by choosing a small positive e in 
the range < e < t + T — t'). Similarly 

A T n(r, 7-0 = A t u;(t) ,g(rO = Q V(r) <?(r0 

ff(r0 / * dy i>(y) k(y) + QA g( T ) 9 {t') , 



k T H^{T,T')=Ag[T)^{T') 

A T H 1 p g ( r ,T') = Ag(T)g(T') ; 
A T H gg (T, T >) =0, 



(53) 
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so that the periodicity equation (|40)) takes the form 

A T G(T,T')=g(T)^(T') (aA + 1^ 

+g(T)g(T')(- + !3A--J dy^(y)k(y)) 
= 0, (54) 

whence 



T.+T 

= P = + ~Kq j ' d v^(y)Hy)- (55) 



Finally, the periodicity requirement does not impose 
any restriction on the coefficient 7, but it immediately 
follows from the gauge fixing condition (p6| . 



T.+T / T„+T 

7 = ^ J dyip(y) k(y) - -^-^ \ j dy < -(//') /,■(//) 

r»+T 

-7^2 y d y uj (y) k (y)- 



(56) 



V. VARIATION OF THE DETERMINANT 



The calculation of the variational term ([37]) is based on 
integration by parts and a systematic use of the Wron- 
skian relation (l20l) together with equations for g and tp. 
Using the obtained Green's function (|42|) we see that var- 
ious terms of (|37j) transform as follows. 

The first term contributes zero because the coincidence 
limit G f (t, t) is vanishing, 

Tr (SF G F {r,r'))=j dr sf G F (t, r') \ t , =t 



(57) 



The contribution of the rest of the terms of (|42|) repre- 
sent integrals over the period T of non-periodic functions, 
which explicitly depend on the choice of the initial inte- 
gration point in 



= <p dr S[ 1 ) G f (t,t) =0. 



dr(...) 



z+T 



dr (...). 



(58) 



To check that in the final answer this dependence can- 
cels out, as it should for the full periodic function 
5F Gf{t,t), we keep z arbitrarily fixed in what follows. 



z+T 



Thus, the contribution of the f2-term transforms as 
Tr^FO^r')) =J dr (uj 6F* g + g SF w) 

= - j dr (ujF~ 5g + gF~ Su?j 

= -j z + dT (ufsg + g^6 g 
+ (w[LJ,5g}+W[g,Su]\ * 

[Z+T 

= drkSg + A T W[u),8g\{z) 

+ A T W[g,S g u]{z). (59) 

where we used the equations Fg — and f"uj = —k 
and their variational corollaries (<5i^) g — — ~F 5g and 

(<5i^)w = —F'SgUJ — 5 g k = —F^SgU. Similarly, the con- 
tributions of H^, H^ g and H gg terms read 

Tr((5F) H^(r,r')) = J dr^sfyj 

dr ip^Sip = A T W[i>,8ip](z) , (60) 



z+T 



Tr((5F H^ g (r, r')) = J + dr (V sf g + g sf ^) 
= - J dr (i>F* 8g + gF* Stp^j 

= & T {W[il>,6g](z) + W\gM]{*i) > (61) 
Tr(6F H gg (T,r')) =J + drgSfg 

= -[ + dTgfSg = A T W[g,Sg}(z) = 0. (62) 

Therefore, the variation (I3T[) takes the following form in 
terms of the monodromies of the set of Wronskians 

z+T 

Tr (SF G(t, t')) =^ JdrkSg 

Z 

+i [A t W[uj, 6g](z) + A T W[g, Su;}(z)) 

+/? (A T W[^, Sg](z) + A T W[g, <^](z)) . (63) 



These monodromies in their turn express via the mon- 
odromy parameter A of the function ip( T ) an d its vari- 
ation. First, in view of the relations W[g,g\ = and 
W[ip,g} = — 1 we have 

A T W , S^} (z) = W [ (V + A g) , 8 + A g) } (z) 
-W[if>,6il>](z) 
= -SA + A W[ip, Sg] + A W[g, Si/j] 
+A 2 W[g,6g]. (64) 
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Second, using the relation W[g, dip] = 5W[g,ip] — 
W[Sg, ip] we have the chain of identities for the 

k T (w[i>,5g] + W[g,8^])(z) 

= A T (W[i/>, Sg] + SW[g, V] - W[Sg, (z) 

= 2W[(1> + Ag),5g](z) - 2W[il>,5g](z) 

= 2AW[g,Sg](z), (65) 

where we used the fact that W[g, tp] = 1 and <5W[<?, ip] = 
0. Similarly 

k T (w[u,8g] + W[g,8 g <4){z) 

= A T (W[u, Sg] + SW[g, w] - W[Sg, w]) (z) 
= 2W[A T u,5g]{z) + 5A T Q(z), 

because W[w, <7](r) = Q(r) and A T SQ = SA T Q. Finally, 
we use the monodromy of uj given by Eq. (|51[) and the fact 
that the monodromy of the function Q(t) defined by f4"5)l 
equals A T Q(z) = Q where Q is defined by Eg. dm) . Thus 
we have for this pair of monodromies 

A r (W[u, Sg] + W[g, Su]) (z) = 2Q W[^(t), Sg](z) 

T,+T 



2{AQ 



dy^(y)k(y))W[gJg](z) + SQ. (66) 



Using in (|63l) the expressions (|55|) for a and (3 and 
the above set of Wronskian monodromies (|M)) - (I56^) we 
eventually find 

Tr(,5F G(r, r')) =2^ + ^ + W{i/;,Sg}(z) 



-W[g,H]{z) 

-- 2 5 In Q + 5 In A . 



(67) 



Here the z-dependent terms proportional to W[g, $g]{z) 
and in W[i>,Sg](z) - W[g,Sip](z) = SW[tfj 7 g] = duti- 
fully cancel out. Therefore, the total variation (l33l) be- 
comes independent of z and gauge independent (as it 
should, c.f. Eq.(|55jl) 



51n(Det*F) = 2 5 In + S In A, (68) 
and we finally get the explicit answer for Det* F 

DeUF = C{T)xA<j)dTg 2 (T). (69) 

This result is valid up to an overall coefficient C(T) 
which is functionally independent of g(r), but can de- 
pend on the period T - the only remaining free param- 
eter. This coefficient function can be determined by the 
zeta-function method for a particular case of the con- 
stant function g(r) = c corresponding to the operator 



F = —d 2 /d,T 2 with the explicit spectrum of eigenfunc- 
tions and respective eigenvalues 

ipo = 1, A = 0, 

<Piu(t) = sin (27rnr/T), (^(t) = cos (2imr/T) , 
A„ = (2vm/T) 2 , n = l,2,.... (70) 

The logarithm of the corresponding restricted determi- 
nant - the product of all nonvanishing eigenvalues regu- 
larized by zeta-function method - equals 



lnDct * = 2 E ln hr 



4 In 



)Cfl(0)-4^(0) = 21nT. (71) 



Here Cr( s ) = Y^=\ n s ^ s the Riemann zeta function 
having the following particular value Cfl(0) — —\ and 
the value of its derivative Cr(0) = — \ hi 27r. On the other 
hand, the basis functions and the monodromy parameter 
A for this operator read 



g(r)=c, jdTg 2 = c 2 T, 
iKt) = -(t-t„), A = -j. 



(72) 



Therefore according to © Det* (-d 2 /dr 2 ) = C*(T)T 2 , 
and the comparison with (|71[) gives the T-independent 
result C(T) = 1, which confirms the unit normalization 
coefficient in (l24l). 



VI. THE MONODROMY ALGORITHM 

From the definition of the monodromy parameter ([23 
it follows that 



#r + T) - t/>(t) 
5(r) 



(73) 



By tending the point r to the one of the roots of <?(r), 
t — ?• tq , r + T — » T2fe , we have 



A = 



1p{T2k) ~ 1p(T ) 




where we took into account that g(r2k) — g( T o)- Thus 
we immediately arrive at the algorithm (|26p -([27 |) . It is 
important that unlike in the construction of the func- 
tion i/j(t) which has to be smooth on S 1 at all roots Tj 
except To (the property that was attained above by a 
special choice of the auxiliary points r*), the derivatives 
of neighboring functions ^(r) in ([2"B)) - ([2"7]) should not 
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necessarily be matched at these junction points. This is 
because the partial contributions Aj to the overall mon- 
odromy parameter A are individually independent of r* , 



dA 

dr* 



0, 



(75) 



which can be easily verified by using a simple relation 
d&i(T)/dT* = -g(T)/g 2 (T*). Thus, the monodromy is 
uniquely defined and independent of the choice of the 
auxiliary points t* necessarily entering the definition of 
functions S^(r) in Eq. (j25|) . 

To derive the manageable integral expression (j2"9")l for 
<?(r) at the roots t± of a zero mode, we note that its 
derivative 



9(t) = g(r) 



1 



dy _^ 

9 2 (y) 9{r) 



(76) 



has two terms which separately diverge when r — > t±. 
To circumvent this difficulty one can add and subtract in 
the integrand of the first term the total derivative term in 
the integration variable y, (l/g(T))(d/dy)(l/g(y)). This 
yields the representation for the time derivative, 



<?(r) = / dy 



9(r) - g(y) 

g 2 (y) 



(77) 



vacuum Casimir energy. Here rap = (37r/4G) 1//2 and A 
are the renormalized Planck mass and cosmological con- 
stant, a' = a/N and the free energy F(rf) is the function 
of 77 - the inverse of the effective temperature, which is 
given by the full period of the conformal time, see ([5]). 
This is a typical boson or fermion sum over field oscilla- 
tors with energies u on a unit 3-sphere. 

Semiclassically the integral © is dominated by the 
saddle points — cosmological instanton solutions of the 
effective Friedmann equation of the Euclidean gravity 
theory 



a' z 1 
« + ^-B 



2a 4 



A 
3~ 



C_ 



c = 



B ] 1 dFjrf) 
2 rrip drj 



(80) 
(81) 



This Friedmann equation is modified by the quantum B- 
term and contains the radiation term C/a , the constant 
C characterizing the sum of the renormalized Casimir 
energy B/2 and the energy of the gas of thermally excited 
particles dF(r))/dr] = XL u>/(e u,ri T 1). 

As shown in [5], these solutions have a fc-fold nature 
(with any integer fc = l,2,3,...). This means that during 
the full period of its Euclidean time T the scale factor 
oscillates k times back and forth between its maximum 
and minimum values a± = a(r±), a_ < a(r) < a+. 



which now allows one to take the limit t = t± directly in 
the integrand, integration remaining convergent because 
the value of the numerator function g{T±)—g{y) and its y- 
derivative, cx g{y), are vanishing at y = t± (cf. Eg. (1181) ). 
This gives ([29]). 



A. CFT driven cosmology: single node case vs 
multiple nodes 

CFT driven cosmology considered in 0, 0, Q is char- 
acterized by the Einstein-Hilbert action with matter 
sources in the form of many field species conformally cou- 
pled to metric. For a large number of these species the 
path integral (j4]) in this model with a good accuracy gen- 
erates the effective action ([5]) with the local Lagrangian 
of minisuperspace variables and the free energy of matter 
fields 



£(a, a!) 



A 3 
- a+ 3 fl 

, a' 2 a' 4 1 
-B I 1 

a 6a 2a 



F(T,)=±^ln(l T e-^. 



(78) 
(79) 



The Lagrangian C(a,a') includes the minisuperspace 
Einstein-Hilbert term, the contributions of the overall 
conformal anomaly of these fields (B > is the coef- 
ficient of its topological Gauss-Bonnet term) and their 



21^^ 



4AC/3), 



(82) 



and forms a kind of a garland of S 1 x S 3 topology with the 
oscillating S 3 section. Therefore, the inverse temperature 
77 - the instanton period in units of the conformal time - 
is given by the integral 



T() + fcT T + 

dr — = 2k dr — 
a J a 

T T_ 



(83) 



over the full period T which is the 2fc-multiple of the 
integral between the two neighboring turning points t± 
of the scale factor history or fc-multiple of one period T 
of such an oscillation, 

d(r±) = 0, T = 2(t+-t_), T = kT. (84) 

Thus, for k > 1 the period T of the periodic functions on 
the functional space of which we calculate the determi- 
nant Det* F is the fc-multiple of the actual period T of 
the zero mode, so that the total period of the instanton 
is composed by glueing together fc such fundamental pe- 
riods. Each of these periods contains two turning points 
of a(r), and since these points coincide with the roots of 
the zero mode g cx d this fc-fold instanton gives rise to fc 
nodes of g(r). 

This situation is a particular case of our general setting 
of Sect. 2, which corresponds to the "equidistant" set of 
roots of the zero mode function © associated with the 
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turning points of the instanton solution. One can check 
that all the properties p^|) -([T^ | of the function g(r) de- 
fined by Eq.© for the La grangian (|78|) are satisfied 6 . 
Moreover, any solution a(r) of this equation is symmetric 
under the reflection of the time variable relative to any of 
its turning points, a(r± + r) = a(r± — r), so that a and 
g oc d are antisymmetric under this reflection. There- 
fore, one root of g can be placed right in the middle of 
the fundamental period of length T, so that all segments 
of time [Tj_i , Tj] connecting the pairs of neighboring roots 
are of the same length and g(r) is antisymmetric under 
reflections with respect to each of these roots. 

For the simplest single- node case with k — 1 which was 
considered in [j| [H, Ho| the location of the roots can be 
parameterized as 



TO = -7U 



7-1=0, 



T 2 = T 4 



(85) 



with the points r and r 2 identified (in order to make the 
time range, — t + < r < r + , to form the circle just as in 
([I])])). Then g(r) is an odd function of r, 



(86) 



and has two first degree zeros at antipodal points of 
this circle r = r_ = and r = ±r + which mark the 
boundaries of two half-periods of the total time range, 
-T/2 < t < T/2. Correspondingly the functions %{t) 
defined by Eq. ([2l))) can be chosen in the form 



W 2 (t) =Wr(r)=g(r) 



dy 
9 2 (y)' 



< t < t+, (87) 



I g 2 (y) * { ~ 



■r), -ti < t < 0, 



with some r* > in the half period range of r. With this 
choice the monodromy constituents Aj become 

Ai = A 2 = -(<£ f «£ f -#_<?•_), (89) 

*±=9(t±), &±=&(t ± ), (90) 

and the total monodromy coefficient reproduces the re- 
sult obtained for a single-node case in [18| (denoted there 
by I) 



k=l 



2Ai = -2(#+lP+ (91) 



Due to the "equidistant" property of the set of roots 
discussed above all terms of the total monodromy param- 
eter coincide with ([8"9"]l. Aj = A/2, and the mon- 
odromy parameter for a multi-node background equals 



A 



fc>i 



= fcA, 



(92) 



The last property (1 1 8 P is guaranteed by the fact that for any 
solution of Eg . (ISOfl all odd order derivatives of a,(r) are vanishing 
at its turning points. 



where A can be calculated by means of Eq. (|9Tj) with the 
function &(t) prescribed on the single half-period of the 
zero mode by Eq. ([87| . 7 



VII. CONCLUSIONS 

Thus we obtained the restricted functional determi- 
nant of a special second order operator (JXJ) , subject to 
periodic boundary conditions with a periodic zero mode 
g = <?(r). This operator F determines the one- loop 
statistical sum for the microcanonical ensemble in cos- 
mology generated by a conformal field theory (CFT). 
This ensemble realizes the concept of cosmological ini- 
tial conditions by generalizing the notion of the no- 
boundary wavefunction of the Universe to the level of 
a special quasi-thermal state which is dominated by in- 
stantons with an oscillating scale factor of their Euclidean 
Friedmann-Robertson- Walker metric. These oscillations 
result in the multi-node nature of the zero mode g(r), 
©J of F, which is gauged out in the path integral for 
the statistical sum by the method of the Faddeev-Popov 
gauge fixing procedure. This effectively leads to the re- 
stricted functional determinant (with the zero eigenvalue 
omitted) . 

In contrast to the previous work [18] we systematically 
used the monodromy method for the calculation of this 
determinant, which allowed us to generalize the previ- 
ously known result for the single-node case [l8[ to the 
general case of multiple nodes (roots) of g(r) within the 
instanton period of the Euclidean time r. The functional 
determinant of F expresses in terms of the monodromy of 
its basis function (|24p . which is obtained in quadratures 
as a sum of contributions of time segments connecting 
neighboring pairs of the zero mode roots within the pe- 
riod range ([26|) - ([27| . The advantage of this representa- 
tion is that the solutions of homogeneous equation explic- 
itly known in quadratures for these different segments 
([231) should not necessarily match with one another to 
form a smooth basis function on the entire period, which 
leaves integration constants r* in (1231) arbitrary. This 
disentangles these segments and makes (|26p an additive 
sum of independent contributions. 

For the single-node case of the CFT driven cosmology 
(k = 1 number of oscillations of the background scale 
factor or Ik = 2 roots of the zero mode) we reproduce 
the result of [llj]. Moreover, in view of the equidistant 



7 One should not think that the monodromy for the fc-node case 
is numerically k times larger than its value for the single-node 
one, because the histories a(r) for multi-fold instantons are not 
obtained by simply matching together k copies of the single-fold 
instanton solution. Simple additivity is violated in view of the 
"bootstrap" equation 1811 relating the amount of radiation C 
on the instanton background to its geometry. For this reason, 
in particular, the on-shell value of the instanton action is not a 
simple multiple of k, but rather tends to zero for k — > oo ||. 



11 



nature of the set of these roots for multi-node case, the 
monodromy becomes a 2fc-multiple of the contribution 
of the half-period of the zero mode g{r) (|89 |) -([92 ]) . Ap- 
plication of this result for a single-fold instanton in [20 ] 
showed smallness of its one-loop part in the statistical 
sum of the CFT driven cosmology. However, preliminary 
analysis for multiple-fold instantons shows the possibil- 
ity of a phase transition, strong coupling domain and also 
peculiar behavior for k — > oo. This limit corresponds to 
the upper bound of the admissible cosmological constant 
range A max = 3/2B - a new quantum gravity scale which 
is of special interest in the microcanonical ensemble of 
the CFT driven cosmology d, ||. Cosmological multi- 
fold instantons of @ are analogous to multi-instanton 
solutions in Yang-Mills theory, but in contrast to the 
Yang-Mills case no dilute gas approximation applies here, 
and their action is not a multiple of the single-fold one 



and, moreover, tends to zero for k — > oo [5| (see footnote 
above). Thus, their contributions are not exponentially 
suppressed (or infinitely enhanced like in 21] in view of 
the negative value of the classical gravitational action), 
and play important role in the full ensemble. Therefore, 
the above results for k > 1 become indispensable for the 
calculation of their preexponential factors to be consid- 
ered in HH. 
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